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A  new  dissociation  model  for  the  direct  simulation  Monte  Carlo  (DSMC)  method  is  formulated 
using  the  threshold  line  concept.  This  model  considers  in  some  detail  the  coupling  between  the 
vibrational  energy  distribution  of  molecules  in  a  gas  and  the  rate  of  dissociation.  For  a  particular 
reaction,  the  new  threshold  line  model  only  requires  determination  of  a  leading  multiplicative 
constant  and  is  therefore  significantly  less  deterministic  than  existing  DSMC  chemistry  models.  The 
new  model  is  evaluated  in  a  heat  bath  under  equilibrium  and  nonequilibrium  conditions.  It  is  also 
employed  to  calculate  hypersonic  shock  waves  of  nitrogen  and  air.  Comparison  of  the  results 
generated  with  the  new  model  is  made  with  existing  experimental  data  and  with  previous 
computational  results.  The  more  detailed  DSMC  implementation  reveals  properties  of  the  threshold 
line  dissociation  model  not  apparent  in  the  previous  continuum  implementation.  ©  1996  American 
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I.  INTRODUCTION 

In  the  present  study  we  focus  on  the  well-known  mecha¬ 
nism  whereby  the  rate  of  dissociation  is  altered  under  condi¬ 
tions  of  thermal  nonequilibrium  behind  shock  waves.  Spe¬ 
cifically,  the  dissociation  rate  is  reduced  when  the  vibrational 
temperature  is  lower  than  the  translational  value.  This 
mechanism  is  often  called  coupled  vibration  dissociation 
(CVD).  It  is  important  to  capture  this  phenomenon  in  com¬ 
putational  models  applied  to  high-temperature  nonequilib¬ 
rium  flows.  For  example,  in  hypersonic  flow  of  vehicles  re¬ 
entering  the  atmosphere  of  the  earth,  chemiical 
nonequilibrium  affects  the  heat  load  and  pitching  moment 
experienced  by  the  body.  In  addition,  the  vibrational  tem¬ 
perature  is  directly  affected  by  the  CVD  model  and  this 
quantity  is  often  used  to  estimate  radiation  emitted  from  hy¬ 
personic  shock  waves. 

The  two  most  popular  approaches  for  performing  nu¬ 
merical  simulations  of  nonequilibrium  flows  are  the  various 
continuum  based  techniques  of  computational  fluid  dynamics 
(CFD),  and  the  particle-based  direct  simulation  Monte  Carlo 
method  (DSMC).  In  CFD  methods,  the  continuum  equations 
of  fluid  dynamics  are  solved  using  finite  difference  or  finite 
volume  formulations.  In  the  DSMC  technique,  a  direct  simu¬ 
lation  of  the  gasdynamics  is  performed  using  particles  that 
move  in  physical  space  and  undergo  intermolecular  colli¬ 
sions.  Clearly,  with  the  DSMC  approach,  there  is  the  oppor¬ 
tunity  to  include  more  detailed  physics  for  collision  phenom¬ 
ena.  Unfortunately,  this  more  direct  approach  is  more 
numerically  intensive  limiting  application  of  the  DSMC 
technique  to  flows  of  relatively  low  density  (e.g,,  Knudsen 
numbers  larger  than  10”"^). 

Several  CVD  models  have  been  proposed  over  the  years 
for  both  CFD  and  DSMC  simulations.^"^  The  models  most 
frequently  employed  in  large-scale  numerical  simulations  of 
nonequilibrium  flows  are  highly  phenomenological.  The  goal 
of  this  paper  is  to  investigate  use  of  a  more  detailed,  less 
phenomenological  CVD  model  for  the  DSMC  technique. 
The  approach  adopted  is  to  use  the  threshold  line  concept 


employed  by  Macheret  and  Rich"^  to  derive  a  continuum 
CVD  model.  The  development  of  such  a  model  for  the  par¬ 
ticle  DSMC  technique  requires  adaptation  of  the  threshold 
line  method  to  a  detailed  microscopic  kinetic  theory  analysis. 
It  will  be  shown  that  this  analysis  reveals  features  of  the 
threshold  line  concept  not  previously  apparent  through  the 
continuum  analysis  of  Ref.  4. 

The  outline  of  the  paper  is  as  follows.  First,  the  standard 
CVD  models  are  described  for  flow  simulations  based  on 
continuum  and  particle  representations  of  the  fluid.  Then,  the 
continuum  CFD  model  proposed  by  Macheret  and  Rich"^  that 
is  based  on  the  threshold  line  concept  is  introduced.  In  the 
present  study,  based  on  work  of  Macheret  et  a  threshold 
line  model  suitable  for  the  particle  approach  is  considered. 
Basic  properties  of  the  model  are  investigated  under  equilib¬ 
rium  and  nonequilibrium  conditions.  Results  obtained  with 
the  new  model  are  compared  with  experimental  measure¬ 
ments  in  hypersonic  shock  waves  of  nitrogen  and  air. 


II.  CVD  MODELS 

As  mentioned  previously,  a  number  of  models  for  simu¬ 
lating  the  coupling  between  the  distribution  of  vibrational 
energy  and  the  dissociation  rate  have  been  proposed.  Some 
of  the  more  popular  models  are  described  in  Refs.  1-6.  Of 
these,  all  but  that  of  Boyd^  are  based  on  continuum  formu¬ 
lations.  In  the  following,  the  CVD  models  presented  in  Refs. 
2  and  6  are  described  in  detail  due  to  their  relatively  wide¬ 
spread  use  in  continuum  and  particle  simulations  of  nonequi¬ 
librium  flows. 

A.  Continuum  approach 

In  continuum  computational  fluid  dynamics  (CFD),  the 
two-temperature  model  of  Park^  is  widely  used.  In  this 
model,  an  Arrhenius  reaction  rate  coefficient  has  its  tempera¬ 
ture  dependence  governed  by 
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where  is  the  translational  temperature,  is  the  vibra¬ 
tional  temperature,  and  ^  is  a  modeling  parameter.  This 
model  is  basically  empirical  but  has  enjoyed  a  fair  degree  of 
success.  It  is  clearly  phenomenological  in  that  it  retains  equi¬ 
librium  behavior  when  the  translational  and  vibrational  tem¬ 
peratures  are  equal  and  results  in  a  reduction  of  the  dissocia¬ 
tion  rate  when  the  vibrational  temperature  is  less  than  the 
translational  value.  The  value  of  the  parameter  s  is  usually 
set  to  0.5  for  dissociation  of  molecular  nitrogen  and  oxygen. 
As  an  indication  of  the  empiricism  involved  with  the  model, 
however,  some  studies  have  recommended  use  of  a  value  of 
0.7.^  The  basic  failing  of  the  model  is  that  there  is  no  way  to 
determine  an  appropriate  value  for  s  except  through  direct 
comparison  with  experimental  data. 


B.  Particle  approach 

In  the  particle-based  direct  simulation  Monte  Carlo 
(DSMC)  method,  the  vibrationally  favored  dissociation 
(VFD)  model  proposed  in  Refs.  6  and  8  is  considered.  In  this 
model,  the  probability  that  a  particular  collision  of  two  par¬ 
ticles  leads  to  dissociation  is  given  by 

/>,=/(£,)  XE^  (2) 

where  /(E^)  is  a  function  of  the  total  collision  energy,  ; 

is  the  vibrational  energy;  and  is  a  modeling  parameter. 
The  total  collision  energy  is  the  sum  of  the  translational, 
rotational,  and  vibrational  energies  of  the  colliding  particles. 
In  Refs.  6  and  8,  it  has  been  determined  empirically  through 
comparison  with  experimental  data  that  appropriate  values 
for  the  parameter  ^  are  2  for  nitrogen,  and  0.5  for  oxygen. 
The  mathematical  form  for  this  model  is  based  on  analytical 
convenience  rather  than  physical  behavior.  By  including  the 
additional  term,  in  which  the  vibrational  energy  is  raised  to 
the  power  coupling  between  vibrational  energy  and  disso¬ 
ciation  rate  is  included  at  a  simple  level.  Particles  with 
higher  vibrational  energies  produce  higher  probabilities  of 
dissociation  than  lower  lying  particles.  At  the  macroscopic 
level,  the  model  produces  the  desired  CVD  behavior,  in  that 
when  the  vibrational  temperature  is  less  than  the  translational 
temperature,  the  dissociation  rate  is  reduced  below  its  equi¬ 
librium  value. 

As  with  the  Park  model  for  continuum  simulations,  the 
disadvantages  of  the  VFD  model  are  its  high  degree  of  em¬ 
piricism  and  the  fact  that  there  is  no  systematic  way  to  de¬ 
termine  (j)  for  each  new  species  of  interest. 


III.  THRESHOLD  LINE  MODELS 

The  threshold  line  concept  is  based  on  the  premise  that  a 
certain  amount  of  translational  energy  in  a  collision  is  re¬ 
quired  above  the  vibrational  energy  of  the  molecule  to  effect 
the  dissociation  reaction.  This  required  energy  is  the  thresh¬ 
old,  F.  While  the  model  is  developed  from  a  collisional 
standpoint,  emphasis  has  been  placed  in  Refs.  4  and  5  on 
continuum  forms  of  the  model.  In  the  following  sections, 
both  continuum  and  particle  forms  of  the  threshold  line  con¬ 
cept  are  discussed. 


A.  Continuum  approach 

In  the  work  of  Macheret  and  Rich,^  it  is  proposed  that  at 
high  vibrational  energy  levels  the  threshold  energy  is  simply 
the  difference  between  the  activation  energy  of  the  reaction 
and  the  molecular  vibrational  energy.  By  detailed  consider¬ 
ation  of  the  collision  dynamics  involved  in  the  dissociation 
process,  it  is  proposed  that  the  threshold  energy  at  low  vi¬ 
brational  levels  is  larger  than  this  simple  difference. 

Using  these  ideas,  continuum  models  expressed  in  terms 
of  translational,  rotational,  and  vibrational  temperatures  have 
been  proposed  by  Macheret  et  al^  to  model  nonequilibrium 
dissociation  processes.  Current  interest  in  the  threshold  line 
concept  has  arisen  in  part  from  the  successful  application  of 
the  continuum-based  model  of  Ref.  4  by  Boyd  et  al^  In  this 
study,  the  threshold  line  model  gave  the  first  successful  pre¬ 
dictions  of  data  from  a  high-altitude,  hypersonic  flight  ex¬ 
periment. 

The  continuum  models  presented  in  Refs.  4  and  5  are 
based  on  integration  of  collision  probabilities  over  equilib¬ 
rium  energy  distributions.  To  obtain  analytical  results,  it  was 
necessary  in  these  studies  to  employ  approximations  in  this 
integration  process.  Specifically,  the  method  of  steepest  de¬ 
scent  and  a  continuous  form  of  the  vibrational  energy  distri¬ 
bution  were  employed. 

A  key  point  about  the  continuum  threshold  line  models 
is  the  manner  in  which  their  behavior  is  constrained  under 
equilibrium  conditions.  The  approach  adopted  is  to  ensure, 
when  all  the  temperatures  are  equal,  that  the  rate  coefficient 
provided  by  the  model  corresponds  to  the  modified  Arrhen¬ 
ius  form: 

kd=aT^  t\i^{-DlkT),  (3) 

where  a  and  are  constants,  k  is  the  Boltzmann  constant, 
and  D  is  the  activation  energy  of  the  reaction.  Hence,  the 
constants  a  and  7]  appear  directly  in  the  continuum  versions 
of  the  threshold  line  model.  One  subtle  effect  of  imposing 
this  condition  on  the  model  is  that  a  specific  temperature 
dependence  is  being  imposed  on  the  probability  of  reaction. 
This  approach  is  not  applicable  to  the  particle  method  that 
requires  dissociation  probabilities  to  be  described  purely  in 
terms  of  collision  energies  as  described  in  the  following  sec¬ 
tion. 

B.  Particle  approach 

The  purpose  of  this  study  is  to  use  the  threshold  concept 
to  develop  a  nonequilibrium  dissociation  model  that  may  be 
employed  directly  in  the  DSMC  technique.  This  requires  the 
formulation  of  a  dissociation  probability  based  on  the  trans¬ 
lational  collision  energy  and  the  rotational  and  vibrational 
energies  of  the  dissociating  molecule.  This  has  been  per¬ 
formed  in  the  work  of  Macheret  et  al}  However,  the  final 
results  of  that  study  are  described  in  terms  of  continuum 
temperatures  rather  than  collision  energies.  For  simplicity, 
the  following  discussions  will  be  concerned  with  dissociation 
of  molecular  nitrogen  through  collisions  with  other  nitrogen 
molecules  or  atoms. 

As  described  in  detail  in  Ref.  5,  the  threshold  line  prob¬ 
ability  of  dissociation  has  distinct  mathematical  forms  for 
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low  and  high  vibrational  levels.  At  low  vibrational  levels, 
there  is  insufficient  velocity  along  the  line  of  centers  of  the 
two  atoms  that  constitute  the  molecule  to  match  the  transla¬ 
tional  velocity  of  the  colliding  atom.  This  is  the  requirement 
that  must  be  met  to  allow  all  of  the  translational  energy  of 
the  atom  to  be  transferred  into  the  vibrational  energy  of  the 
molecule  and  thus  cause  dissociation.  Hence,  at  low  vibra¬ 
tional  energies,  the  threshold  energy  required  to  cause  disso¬ 
ciation  is  greater  than  the  difference  between  the  activation 
energy  and  the  translational  collision  energy.  As  proposed  in 
Ref.  4,  the  boundary  between  high  and  low  vibrational  levels 
is  determined  using  the  following  parameter: 


a= 


(4) 


where  m  is  the  mass  of  one  of  the  two  atoms  in  the  diatomic 
molecule  that  is  dissociating,  and  M  is  the  mass  of  the  col¬ 
lision  partner.  In  the  case  of  diatom-diatom  collisions,  M  is 
the  mass  of  one  of  the  atoms  in  the  second  molecule,  and 
therefore  a  is  equal  to  0.25  for  both  N2-N2  and  N2-N  col¬ 
lisions. 

In  the  following  expressions,  is  the  translational  col¬ 
lision  energy,  and  and  E^  are  the  rotational  and  vibra¬ 
tional  energies  of  the  reacting  molecule.  In  this  study,  the 
vibrational  energy  distribution  is  modeled  as  a  simple  har¬ 
monic  oscillator  using  the  approach  of  Bergemann  and 
Boyd.^®  The  threshold  line  model  specifies  low  vibrational 
energy  levels  as  those  for  which  where  D*  is  a 

modified  activation  energy  that  will  be  defined  later.  Using 
these  ideas,  it  may  be  shown^^  that  the  dissociation  probabil¬ 
ity  of  low  vibrational  levels  for  diatom-atom  collisions  is 


_ S^l2{E-Ff^ _ 

(5) 


where  A,  is  a  constant  of  proportionality  for  reaction  i  that 
does  not  appear  in  the  continuum  formulation.  It  is  necessary 
here  in  order  to  calibrate  the  DSMC  model  against  specified 
rate  coefficients.  In  the  continuum  model,  A,  has  some  un¬ 
known  functional  form  that  ensures  the  Arrhenius  rate  coef¬ 
ficient  is  recovered  at  equilibrium.  For  low  vibrational  levels, 
the  threshold  energy  is  given  as 


I  —  a 


(6) 


At  high  vibrational  levels  {E^>aD^),  the  threshold  en¬ 
ergy  is  simply 

F=D^-E^.  (7) 

Following  the  procedures  outlined  in  Ref.  5,  the  dissociation 
probability  at  the  high  vibrational  levels  for  diatom-atom 
collisions  is  derived  here  as 

„  .  /l  +  V^  (£ -F)2 

F^  —  Aj'A/  I -  -  I - —5  (8) 

V  1-4^  4'rr^F^'^^lE^-aD* 

where  the  constant  of  proportionality  A  ^  is  the  same  as  that  in 
Eq.  (5). 


The  participation  of  rotational  energy  in  the  reaction  pro¬ 
cess  is  included  by  proposing  an  effective  activation  energy 
of  the  form^ 

2El^^ 

2)*  =  /)-^+——,  (9) 

where  D  is  the  dissociation  energy.  At  very  high  rotational 
energies,  the  activation  energy  is  reduced  due  to  the  strong 
centrifugal  force  induced  on  the  molecule.  This  aspect  of  the 
model  reduces  the  dissociation  rate  in  comparison  with  the 
equilibrium  value  when  the  rotational  temperature  is  less 
than  the  translational  value.  This  effect  is  much  smaller  than 
the  coupling  with  the  vibrational  energy  mode. 

In  Ref.  5,  separate  threshold  line  models  are  presented 
for  diatom— diatom  and  diatom- atom  collisions.  In  the 
present  study,  however,  the  above  probability  model  (strictly 
valid  only  for  diatom- atom  collisions)  is  employed  for  both 
^^2-^2  N2-N  dissociation  reactions.  This  approach  is 
adopted  due  to  the  requirement  for  a  relatively  simple  math¬ 
ematical  expression  for  use  in  the  numerically  intensive 
DSMC  technique. 

As  with  the  continuum  model,  the  particle  based  thresh¬ 
old  line  probability  must  at  the  macroscopic  level  recover  a 
specified  Arrhenius  rate  coefficient  under  equilibrium  condi¬ 
tions.  Unlike  the  continuum  model,  however,  the  Arrhenius 
constants  a  and  77  do  not  appear  in  the  model.  Instead,  the 
constant  of  proportionality,  A^-,  must  be  obtained  through 
calibration  of  the  model  against  rate  data  by  performing  test 
simulations  under  equilibrium  conditions.  Hence,  with  this 
approach,  it  is  not  possible  to  guarantee  that  the  model  will 
produce  a  specific  temperature  dependence  for  the  dissocia¬ 
tion  rate.  This  could  perhaps  be  achieved  through  inclusion 
of  further  dependencies  of  the  dissociation  probability  on  the 
various  energies  involved  in  the  collision.  This  lies  beyond 
the  scope  of  the  present  study. 

IV.  RESULTS 

The  properties  of  the  threshold  line  dissociation  model 
for  the  DSMC  technique  are  first  assessed  in  simple  heat 
bath  conditions.  These  are  used  to  calibrate  the  model 
against  experimental  measurements  of  dissociation  rate  coef¬ 
ficients.  In  addition,  comparison  of  the  DSMC  model  with 
the  continuum  threshold  line  model  is  conducted.  The  cali¬ 
brated  DSMC  models  are  then  applied  to  normal  shock 
waves  of  nitrogen,  and  to  a  bow  shock  wave  in  air.  In  the 
latter  cases,  direct  comparison  of  simulated  results  with 
available  experimental  data  is  undertaken. 

A.  Equilibrium  heat  bath 

The  DSMC  threshold  line  model  is  used  to  calculate 
dissociation  rate  coefficients  under  conditions  where  the 
translational,  rotational,  and  vibrational  modes  are  in  thermal 
equilibrium  with  Boltzmann  distributions.  An  isothermal 
heat  bath  is  simulated  that  consists  of  50  000  particles.  En¬ 
ergy  is  exchanged  between  the  various  modes  during  colli¬ 
sion,  but  chemical  reactions  are  not  processed.  Instead,  the 
average  dissociation  probability  is  evaluated  over  all  colli¬ 
sions  and  then  converted  into  a  rate  coefficient.  Results  ob- 


Phys.  Fluids,  Vol.  8,  No.  5,  May  1996 


lain  D.  Boyd  1295 


FIG.  1.  Dissociation  rate  coefficient  for  N2-N2  under  conditions  of  thermal 
equilibrium. 


tained  in  this  manner  are  compared  in  Fig.  1  over  a  wide 
temperature  range  with  the  following  rate  coefficient  of 
Byron^^  for  the  N2-N2  reaction; 

exp(-  113  200/7). 

The  units  are  cm^/mol/s.  As  discussed  by  Boyd,^  the  data  of 
Byron  are  employed  as  they  were  measured  under  conditions 
of  thermal  equilibrium.  It  is  quite  remarkable  that  the  DSMC 
model  provides  a  temperature  dependence  that  is  so  close  to 
the  Byron  rate.  It  should  be  noted  that  the  leading  constant  in 
the  DSMC  model  is  explicitly  chosen  to  give  agreement  with 
the  Byron  rate  at  7=  10  000  K.  In  this  case,  =  5. 

In  Fig.  2  the  dissociation  probabilities  predicted  by  the 
model  for  N2-N2  collisions  at  several  different  vibrational 
energy  levels  are  plotted  as  a  function  of  translational  colli¬ 
sion  energy  (expressed  as  a  temperature  by  dividing  by  2 A:). 
It  is  interesting  to  note  that  there  is  a  region  for  all  vibra¬ 
tional  levels  where  the  dissociation  probability  exceeds 
unity.  In  practice,  when  using  a  probability  model  in  a 
DSMC  calculation,  the  reaction  rate  simulated  will  employ  a 
maximum  reaction  probability  of  unity.  That  is,  even  if  the 
reaction  probability  is  equal  to  2.0,  only  one  reaction  is  pro¬ 
cessed.  Thus,  the  dissociation  rates  shown  in  Fig.  1  were 
generated  by  limiting  the  dissociation  probability  to  a  maxi¬ 
mum  of  unity. 
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FIG.  3.  Dissociation  rate  coefficient  for  N2-N  under  conditions  of  thermal 
equilibrium. 


The  behavior  of  the  DSMC  model  for  N2~-N  collisions  is 
shown  in  Fig.  3.  In  this  case,  the  selection  of  rate  data  is 
more  difficult.  The  rate  coefficients  considered  for  this  reac¬ 
tion  are  those  reported  by  Byron  as 

'*N2-n=7.17X  exp(-  113  200/7), 

by  Cary^^  as 

rN^_N=2.7lX  10"^7‘‘-®  exp(-  113  200/7), 
and  by  Ragan  and  Appleton^"^  as 

'•n2-n=  2.76X  10"^7"  exp(  - 1 13  200/7) . 

In  all  cases  the  units  are  cmVmol/s.  The  DSMC  results  are 
obtained  using  a  leading  constant  in  the  threshold  line  model 
of  “  14.  Good  agreement  is  found  with  Refs.  13  and 
14  at  low  temperature,  and  better  agreement  with  Ref.  12  at 
high  temperature. 

Results  for  O2-O2  and  O2-O  dissociation  reactions  are 
shown  in  Figs.  4  and  5,  respectively.  The  Arrhenius  rates  in 
each  case  are  those  reported  by  Byron: 

ro2-02  =  3.27X  10"^7"‘-^  exp(-59  400/7) 


FIG.  2.  Probability  of  dissociation  at  various  vibrational  energy  levels  as  a 
function  of  translational  collision  energy  for  the  N2-N2  reaction. 


FIG.  4.  Dissociation  rate  coefficient  for  O2-O2  under  conditions  of  thermal 
equilibrium. 
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FIG.  5.  Dissociation  rate  coefficient  for  O2-O  under  conditions  of  thermal 
equilibrium. 

and 

'■02-0=  1-95X  exp(-59  400/7’), 

where  the  units  are  cm^/mol/s.  The  threshold  line  model  is 
calibrated  for  each  reaction  at  a  temperature  of  10  000  K, 
giving  leading  constants  of  ^02-02  “  and  Ao^-o  =  ^2. 
Consistent  with  the  results  obtained  for  nitrogen  dissociation, 
the  DSMC  threshold  line  model  for  oxygen  gives  better 
agreement  with  an  Arrhenius  rate  coefficient  having  a  tem¬ 
perature  exponent  of  —0.5. 

The  key  aspect  of  all  CVD  models  is  the  manner  in 
which  molecules  are  selected  preferentially  for  dissociation 
based  on  their  vibrational  energy.  Further  equilibrium  simu¬ 
lations  are  performed  to  analyze  how  this  occurs  at  the  col¬ 
lision  level  in  the  threshold  line  model.  Once  again,  several 
million  collisions  are  simulated  using  100  000  particles  in 
equilibrium  at  10  000  K.  Using  an  accept-reject  approach, 
individual  collisions  that  give  rise  to  dissociation  are  identi¬ 
fied.  For  each  of  these  reactive  collisions,  the  vibrational 
energy  of  the  reacting  molecule  is  recorded.  The  vibrational 
energy  distributions  of  reactive  molecules  obtained  in  this 
manner  with  the  threshold  line  and  VFD  dissociation  models 
for  N2-N2  collisions  are  compared  in  Fig.  6.  The  distribu¬ 
tions  produced  by  the  two  models  are  very  different.  By 


FIG.  6.  Distribution  of  vibrational  energy  for  reacting  nitrogen  molecules 
under  equilibrium  conditions  at  10  000  K. 


FIG.  7,  Dissociation  rate  coefficient  for  N2-N2  as  a  function  of  vibrational 
temperature  at  7,  =  7^=  10  000  K. 

comparison  with  the  VFD  model,  the  threshold  line  model 
dissociates  many  more  molecules  at  low  vibrational  levels. 
This  is  a  significant  effect  that  will  become  important  in 
nonequilibrium  flow  conditions.  Note  that  the  border  be¬ 
tween  low  and  high  vibrational  levels  in  the  threshold  line 
model  occurs  for  N2  at  a  vibrational  quantum  number  of  8. 
At  this  point,  the  vibrational  energy  distribution  of  reacting 
molecules  undergoes  a  distinct  but  smooth  change  in  profile. 
While  the  minimum  in  the  vibrational  energy  distribution  of 
reacting  particles  shown  in  Fig.  6  is  an  interesting  property 
of  the  model,  it  is  not,  unfortunately,  feasible  to  measure  this 
energy  distribution  experimentally. 

B.  Nonequilibrium  heat  bath 

Comparison  is  next  made  of  the  DSMC  threshold  line 
model  with  the  continuum  model  of  Macheret  et  al}  under 
nonequilibrium  conditions.  For  simplicity,  the  translational 
and  rotational  modes  are  equilibrated.  Only  the  effect  of  vi¬ 
brational  nonequilibrium  is  considered.  Once  again,  heat 
bath  simulations  are  performed.  In  this  case,  the  translational 
and  rotational  modes  are  given  a  temperature  of  10  000  K. 
For  each  of  the  vibrational  levels  of  the  harmonic  oscillator 
that  lie  below  the  dissociation  energy,  a  large  number  of 
collisions  are  calculated  to  determine  the  state-dependent 
dissociation  rate.  These  rates  are  then  summed  over  the  Bolt¬ 
zmann  populations  at  different  vibrational  temperatures.  The 
results  of  this  analysis  are  shown  in  Fig.  7  for  the  N2-N2 
reaction.  Comparison  is  made  with  the  continuum  threshold 
line  model.  For  reference,  the  dissociation  rate  produced  by 
the  Park  model  using  5  =  0.5  is  also  shown.  This  illustrates 
the  key  aspect  of  the  threshold  line  models,  which  is  that 
they  predict  substantially  higher  dissociation  rates  under  con¬ 
ditions  of  strong  thermal  nonequilibrium. 

Concentrating  on  the  two  threshold  line  profiles  in  Fig. 
7,  it  is  seen  that  they  give  good  agreement  at  the  lower  and 
upper  values  of  vibrational  temperature.  However,  the 
DSMC  model  predicts  a  higher  dissociation  rate  for  the  in¬ 
termediate  values.  The  reason  for  this  difference  is  partly 
explained  through  consideration  of  Fig.  8.  This  shows  the 
contributions  to  the  overall  dissociation  rate  of  the  lower 
{E^<aD)  and  higher  {E^>aD)  vibrational  levels.  The 
DSMC  model  generally  shows  higher  dissociation  rates  from 
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FIG.  8.  Contributions  of  high  and  low  vibrational  levels  to  the  dissociation 
rate  coefficient  for  N2-N2  as  a  function  of  vibrational  temperature  at 
T,=  T  =10  000  K. 

both  low  and  high  vibrational  levels.  It  is  assumed  that  these 
differences  arise  from  the  approximations  made  in  obtaining 
analytical  expressions  for  the  continuum  model.  In  particu¬ 
lar,  integrals  are  evaluated  using  the  method  of  steepest  de¬ 
scent,  which  can  give  rise  to  significant  inaccuracies.  In  ad¬ 
dition,  the  continuum  model  assumes  a  continuous 
distribution  of  vibrational  energy  whereas  the  DSMC  ap¬ 
proach  employs  the  truncated  harmonic  oscillator. 

C.  Normal  shock  waves  of  nitrogen 

In  this  section,  the  DSMC  threshold  model  is  evaluated 
through  direct  comparison  of  results  obtained  with  the 
threshold  line  models  and  existing  experimental  data  for  hy¬ 
personic  shock  waves  of  nitrogen.  The  data  is  taken  from 
Kewley  and  Homung,^^  who  report  density  profiles  obtained 
from  interferograms.  Two  cases  with  different  upstream  ve¬ 
locities  are  considered.  The  DSMC  simulations  of  these  one¬ 
dimensional  shocks  employ  several  thousand  cells  and  about 
100  000  particles.  In  the  first  case,  the  upstream  velocity  and 
density  are  4.8  km/s  and  4.67X10“^  kg/m^,  respectively. 
These  conditions  give  a  post- shock  mole  fraction  of  atomic 
nitrogen  of  just  0.15.  Thus,  for  this  case,  the  density  profile  is 
primarily  characterized  by  N2-N2  dissociation  reactions.  The 
DSMC  results  computed  with  the  threshold  line  model  com¬ 
pare  very  favorably  with  the  experimental  data,  as  indicated 
in  Fig.  9.  The  second  case  considers  an  upstream  velocity 
and  density  of  7.3  km/s  and  7.48X10"^  kg/m^,  respectively. 
Under  these  more  energetic  conditions,  the  post-shock  mole 
fraction  of  atomic  nitrogen  is  0.60.  Thus,  in  this  case,  the 
density  rise  across  the  shock  is  dominated  by  N2-N  disso¬ 
ciation  reactions.  The  profile  obtained  with  the  threshold  line 
model  is  compared  with  the  experimental  data  in  Fig.  10. 
Once  again,  good  agreement  is  obtained  between  simulation 
and  measurement.  The  level  of  agreement  obtained  with  the 
threshold  line  models  and  the  experimental  data  shown  in 
Figs.  9  and  10  is  nearly  identical  to  that  reported  in  Ref.  6  for 
the  VFD  chemistry  model.  It  should  be  noted,  however,  that 
the  agreement  for  the  VFD  model  was  only  obtained  through 
appropriate  selection  of  the  phenomenological  parameter  cf) 
in  Eq.  (2). 


FIG.  9.  Ratio  of  local  to  upstream  density  in  a  hypersonic  shock  wave  of 
nitrogen:  t/j  =  4.8  km/s. 

D.  Bow  shock  wave  of  air 

To  evaluate  the  threshold  line  models  for  oxygen  disso¬ 
ciation,  the  BSUV-2  flight  experiment  is  considered.  As 
mentioned  earlier,  favorable  comparisons  of  results  obtained 
with  the  continuum  threshold  line  models  and  experimental 
data  from  this  flight  were  reported  in  Ref.  9.  The  solutions  at 
the  highest  altitudes  considered  at  85  and  88  km  were  ob¬ 
tained  in  Ref.  9  using  a  two-step  process.  First,  a  DSMC 
solution  of  the  flow  field  was  generated  that  considered  only 
the  major  species  of  N2  and  O2.  Using  this  as  input,  the 
species  conservation  equations  were  then  solved  for  the  mi¬ 
nor  species  of  NO,  N,  and  O  with  a  finite  volume  overlay 
technique.  Within  this  continuum  overlay  scheme,  produc¬ 
tion  of  atomic  nitrogen  and  oxygen  through  dissociation  re¬ 
actions  was  modeled  using  the  continuum  threshold  line 
model  of  Macheret  and  Rich."^  An  overlay  approach  is  appro¬ 
priate,  as  the  concentrations  of  the  minor  species  are  several 
orders  of  magnitude  lower  than  those  of  the  major  species. 

In  the  present  work,  the  particle  threshold  line  models 
for  dissociation  of  nitrogen  and  oxygen  are  applied  to  the 
flow  field  generated  by  the  BSUV-2  vehicle  at  a  velocity  of 
5.1  km/s  and  an  altitude  of  85  km.  The  geometry  is  approxi¬ 
mated  by  the  front  portion  of  a  sphere  of  radius  10  cm  rep- 


FIG.  10.  Ratio  of  local  to  upstream  density  in  a  hypersonic  shock  wave  of 
nitrogen:  C/i  =  7.3  km/s. 
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FIG.  11.  Temperatures  along  the  stagnation  streamline  of  BSUV-2  at  85  km 
altitude. 


resenting  the  nose  of  the  vehicle.  The  temperature  profiles 
for  the  translational,  rotational,  and  vibrational  modes  along 
the  stagnation  streamline  for  this  case  are  shown  in  Fig.  11. 
Clearly,  a  very  strong  degree  of  thermal  nonequilibrium  ex¬ 
ists.  In  Fig.  12,  concentrations  of  atomic  oxygen  and  nitric 
oxide  along  the  stagnation  streamline  are  shown.  Results  ob¬ 
tained  previously  using  the  continuum  threshold  line  model 
in  the  continuum  overlay  scheme  (MR-CFD)  are  compared 
with  ±ose  computed  using  the  particle  models  presented  ear¬ 
lier  in  this  study  (MR-DSMC).  Excellent  agreement  is  ob¬ 
tained  for  the  two  solutions.  In  the  present  study,  the  DSMC 
computation  of  these  minor  species  is  handled  directly,  with¬ 
out  the  need  for  an  overlay  computation  using  a  weighting 
scheme,  and  this  accounts  for  the  statistical  scatter  seen  par¬ 
ticularly  in  the  NO  results.  The  agreement  obtained  is  sig¬ 
nificant,  as  the  previous  continuum  results  provided  excellent 
agreement  with  experimental  measurements  of  radiative 
emission  of  both  atomic  oxygen  and  nitric  oxide,  as  reported 
in  Ref.  9, 


FIG.  12.  Concentrations  of  atomic  oxygen  and  nitric  oxide  for  the  stagna¬ 
tion  streamline  of  BSUV-2  at  85  km  altitude:  Comparison  of  current  DSMC 
results  with  CFD  overlay  results  of  Ref.  9. 


V.  CONCLUDING  REMARKS 

A  threshold  line  dissociation  model  for  the  DSMC  tech¬ 
nique  has  been  formulated  that  determines  the  probability  of 
reaction  based  on  the  energy  of  each  collision.  For  a  particu¬ 
lar  dissociation  reaction,  the  model  only  requires  determina¬ 
tion  of  a  leading  constant  through  comparison  with  an 
Arrhenius  rate  coefficient.  The  model  was  found  to  repro¬ 
duce  the  temperature  dependence  of  measured  rate  data  un¬ 
der  equilibrium  conditions.  In  agreement  with  the  original 
continuum  threshold  line  model,  the  new  particle-based 
threshold  line  model  predicted  dissociation  rates  under  con¬ 
ditions  of  strong  thermal  nonequilibrium  that  are  orders  of 
magnitude  higher  than  those  predicted  by  existing  dissocia¬ 
tion  coupling  models.  In  application  of  the  new  model  to 
hypersonic  shock  waves  of  nitrogen,  excellent  agreement 
was  obtained  with  experimental  measurements  of  density. 
For  hypersonic,  dissociating  air  at  high-altitude  conditions, 
the  new  models  for  nitrogen  and  oxygen  dissociation  were 
applied  and  found  to  give  excellent  agreement  with  previous 
computational  results  obtained  with  the  continuum  threshold 
line  model. 

An  additional  aspect  of  the  DSMC  implementation  of 
the  threshold  line  concept  for  dissociation  was  that  properties 
of  the  model  were  revealed  that  were  not  apparent  from  the 
previous  continuum  implementation.  First,  a  minimum  in  the 
vibrational  energy  distribution  of  dissociating  molecules  was 
obtained.  Second,  the  rate  of  dissociation  under  thermal  non¬ 
equilibrium  conditions  was  found  to  reach  an  asymptotic 
value  as  the  vibrational  temperature  tended  to  zero. 

These  studies  have  demonstrated  the  physical  accuracy 
of  the  threshold  line  models  presented  in  this  investigation, 
and  their  use  is  recommended  over  previous  phenomenologi¬ 
cal  DSMC  models  for  vibration-dissociation  coupling. 
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